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Abstract 

This  report  is  a  summary  of  a  one-year  effort  on  modeling  air  plasma  actuators  by  the 
Pi’s  research  team  and  collaborators.  Plasma  actuators  have  become  the  key  enabler  for 
boundary  layer  flow  control  especially  when  the  device  needs  to  be  surface  compliant. 
While  the  effects  of  these  actuators  are  striking  for  low  speeds,  their  efficacy  becomes 
very  limited  for  a  wide  range  of  flow  speeds  due  to  the  lack  of  understanding  of  the 
multiphysics  controlling  adjustable  authority  of  such  devices.  In  our  previous  effort  we 
have  demonstrated  the  prediction  capability  of  module-based  multiscale  ionized  gas 
(MIG)  flow  finite-element  code  for  mitigating  inert  gas  flow  separation  using  rf-driven 
dielectric  barrier  discharge.  In  this  effort  we:  (1)  develop  multi-dimensional  first 
principles  based  ^/Chair  chemistry’  models  for  the  non-equilibnum  real  gas  discharge, 
and  (2)  implement  it  in  a  finite  element  framework  to  predict  plasma  discharge 
characteristics  and  its  associated  electric  force  generation  over  a  three-dimensional 
conical  forebody  geometry.  To  our  knowledge,  such  detailed  plasma  kinetics  based  effort 
has  not  been  reported  before.  During  the  development  of  this  project  we  have  worked  in 
close  collaboration  with  RBAC/AFRL  researchers  to  integrate  MIG  results  with  FDL3D1 
code  and  to  simulate  electrodynamic  mitigation  over  practical  surface. 


20090723682 


I.  INTRODUCTION 


The  experiments  and  numerical  predictions  reported  lately  show  the  ability  of 
plasma  actuators  to  impart  body  force  inside  the  boundary  layer  of  a  fluid  in  the  close 
vicinity  of  a  surface.  For  plasma  actuators  with  spatially  displaced  electrodes  the  airflow 
inducement  is  close  to  the  dielectric  surface  with  a  time-averaged  flow  velocity 
measured  between  1-3  m/s.  These  actuators  exhibit  several  benefits  in  active  flow  control 
applications,  including  absence  of  moving  parts,  rapid  on-off  deployment  and  attractive 
self-limiting  characteristics.  Figure  1  represents  a  standard  two-dimensional  geometry 
which  has  been  extensively  investigated  for  boundary  layer  flow  control  at  low  speeds. 
The  asymmetric  positioning  of  the  electrodes  helps  directional  bias  of  the  majority 
electric  force  and  hence  inducement  of  a  wall  jet.  The  induced  force  is  local  (within  a 
few  mm)  and  dissipates  rapidly  outside  the  domain  of  influence. 


Figure  1.  Standard  linear  plasma  actuator- with  primarily  two-dimensional  discharge  and 
induced  velocity. 

The  discharge  may  be  characterized  as  transient  with  a  time  scale  of  the  order  of  a 
fewr  microseconds  while  that  for  the  bulk  gas  flow  is  in  milliseconds.  This  allows  the 
fluid  to  respond  relatively  instantaneously  to  the  electrical  inputs.  Such  force 
successfully  controls  flow  separation  on  airfoils  at  a  high  angle-of-attack,  increases  the 
lift  and/or  reduces  the  drag  of  airfoils  and  fuselages.  Despite  experimental  success, 
applications  of  such  actuators  are  limited  to  low  speeds  due  to  an  incomplete 
understanding  of  the  plasma-gas  coupling  in  all  scales  and  dimensions.  For  example,  to 
our  knowledge,  the  reported  numerical  efforts  are  either  one  or  two-dimensional,  and  the 


experimental  efforts  lack  temporal  and  spatial  resolution  of  data.  Thus  it  is  critical  to 
identify  the  three-dimensional  effects,  if  any,  of  force  generation  due  to  an  operational 
plasma  actuator. 

The  fundamental  principles  of  plasma  actuators  are  continually  being  explored  by 
both  experimental  techniques  and  numerical  modeling.  Recently,  the  dc  corona 
discharges  and  ac  dielectric  barrier  discharges  have  been  investigated  under  a  range  of 
operating  conditions.  In  experimental  studies  the  V-dot  probe  technique  on  actuators  in  a 
large  vacuum  chamber  shows  that  the  applied  potential  across  the  electrodes  is  twice  of 
the  voltage  drop  across  the  plasma  generated  10.  It  has  been  identified  that  the  electric 
force  generation  is  linearly  related  to  the  surrounding  air  pressure  .  Adjustment  of  the 
actuator  geometry,  dielectric  materials,  power  frequency,  and  rms  voltage  show 
improvement  in  power  transferred  to  the  neutral  gas  flow  by  ion-neutral  collisions  l2.  For 
example,  the  electrically  induced  body  force  increases  with  up  to  the  fourth  power  of 
voltage  for  a  constant  frequency  1  The  induced  jet  velocity  thus  improves  with  applied 
voltage  4  \  Experiments  show  that  the  actuator  power  may  be  reduced  by  90  percent 
and  its  durability  can  be  improved  by  operating  it  in  an  unsteady  manner  .  However 
confusion  remains  about  the  influence  of  frequency.  In  one  study  the  body  force  was 
found  less  biased  by  the  input  frequency  n,  while  others  found  it  to  be  proportional  to 
frequency  for  a  constant  voltage  1  The  influence  of  plasma  actuators  are  limited  by 
different  kinds  of  power  losses  such  as  reactive  power  losses  due  to  impedance  mismatch 
between  the  power  supply  and  the  actuator,  dielectric  heating,  plasma  joule  heating  which 
need  to  be  minimized  for  optimum  performance  1  lv 

Specific  modality  of  discharge  whether  it  is  glow,  Townsend,  or  streamers  will 
help  determine  how  the  body  force  is  transported  from  charge  to  neutral  species  19.  The 
effect  of  this  body  force  in  controlling  three  dimensional  laminar  and  transitional  flow 
structures  have  been  studied  utilizing  both  phenomenological  and  first-principles 
plasma  kinetics.  2  Recent  development  of  a  loosely  coupled  fluid  plasma  formulation  22 
has  shown  reasonable  success  in  describing  radio  frequency  (rf)  powered 
electrohydrodynamic  (EHD)  control  of  NACA  wing  stall,  laminar  separation  over  a 
ramp,  boundary  layer  transition  on  a  flat  plate,  and  turbulent  separation  over  a  wall- 
mounted  hump.  Unsteady  flow  actuation  with  a  duty  cycle  seemed  to  perform  better  than 


continuous  operation  of  the  rf  plasma  actuator.  These  results  underscore  the  importance 
of  accurately  resolving  the  unsteady  body  force  and  understanding  its  transition  and 
turbulence  enhancement  mechanisms  for  effective  plasma  actuation. 

Numerical  simulations  of  an  asymmetrically  arranged  dielectric  barrier  discharge 
actuator  employing  multispecies  hydrodynamics  14,23  and  Monte  Carlo  particle-in-cell  and 
direct  simulation  Monte  Carlo  methods  l(  by  independent  research  groups  have  found 
that  ionization  is  not  equal  during  positive  and  negative  part  of  the  alternating  cycle 
producing  a  net  force  in  one  direction  Also,  power  losses  due  to  inadequate  impedance 
matching  of  the  power  supply  to  the  actuator,  dielectric  heating,  and  power  required  to 
maintain  the  atmospheric  pressure  plasma  have  been  reduced  for  better  performance  of 
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the  actuator.  * 

Plasma  actuator  was  also  treated  as  a  heat  and  momentum  input  into  a  small 
control  volume  2.  Reported  measurements  indicate  that  the  momentum  coupling  between 
the  charged  particles  and  the  neutral  particles  occurs  on  timescales  much  shorter  than  that 
for  the  bulk  fluid  motion.  The  description  of  plasma-fluid  interaction  involves  interplay 
between  electron,  ion  and  fluid  time  scales.  The  electron  plasma  frequency  is  -104  Hz, 
ion  plasma  frequency  is  ~106  Hz  and  that  of  the  fluid  (speed/domain)  frequency  is  ~10  ' 
Hz.  As  the  speed  increases,  for  example,  with  Mach  number,  the  flow  frequency 
increases  and  at  very  high  speeds  the  ion  plasma  frequency  may  become  comparable  to 
the  flow'  frequency.  At  hypersonic  flows  non-equilibrium  chemical  kinetics  also  affects 
the  energy  transfer  mechanisms.  Another  mode  of  momentum  coupling  can  be  achieved 
using  wall  temperature  injection  influencing  the  thermal  boundary  layer  through 
electrical  or  other  means.  Theoretical  details  of  such  coupling  will  be  explored  through  an 
experiment-integrated  numerical  study  in  the  near  future. 

An  important  factor  in  understanding  the  force  mechanism  in  plasma  actuators  is 
the  fundamental  influence  of  various  charged  species  in  the  working  gas.  PPs  group  has 
introduced  a  detailed  numerical  model  of  plasma  actuator  with  eight  species  air  chemistry 
with  nitrogen  and  oxygen  mixture  has  found  effect  of  different  charged  species  into  the 
generation  of  clectrodynamic  force  that  controls  flow  of  the  surrounding  fluid 


As  a  specific  application  for  such  model,  a  three-dimensional  study  of  the 
discharge  characteristic  over  a  conical  forebody  geometry  is  chosen.  The  flight  vehicle 
forebody  vortex  symmetry  breaking  and  control  of  resulting  yaw  departure  by  active  and 
passive  means  have  been  studied  by  many  researchers  2t  .  Such  study  includes 
numerical  simulation,  analytical  interpretation  and  experiments.  For  example,  flow 
control  using  flow  injection  and  surface  heating  has  been  investigated  for  asymmetric 
flows  around  circular  cones  by  using  computational  solution  of  the  unsteady, 
compressible  full  Navier-Stokes  equations.  2  Various  computational  aspects  of  two- 
dimensional  modeling  have  already  been  extensively  explored  2(\  Passive  vortical  flow 
control  on  a  wing-body  combination  using  tangential  blowing  has  also  been  studied. 
The  feasibility  of  using  forebody  tangential  blowing  to  control  the  roll-yaw  motion  of  a 
wind-tunnel  model  at  high  angles  of  attack  has  been  showm  experimentally.  *  Analytical 
modeling  for  a  circular  cone  for  determining  the  stability  of  point  vortices  to  small 
symmetric  and  asymmetric  displacements  were  qualitatively  validated  against 
experimental  results  ~9. 

As  an  alternative  to  blowing,  surface  discharge  plasmas  may  be  employed  to 
actively  control  or  eliminate  vortical  asymmetry  about  the  nose  of  a  conical  forcbody  at 
an  angle  of  incidence.  High-frequency  spark  discharge  plasma  was  used  to  control 
asymmetric  vortical  structures  on  the  conical  forebody  30  Results  of  wind  tunnel 
experiments  show'  that  plasma  actuators  distributed  near  How  separation  lines  may  be 
used  for  successful  control  or  elimination  of  vortex  asymmetry.  A  tw'o-dimensional 
theoretical  study  for  flow  control  over  a  conical  section  was  reported  by  the  PTs  group 
using  a  single  barrier  discharge  actuator  and  a  set  of  dc  plasma  actuators  Results  show' 
insufficient  control  of  the  flow  on  the  entire  length  of  the  forebody  due  to  a  single  barrier 
discharge  actuator.  Flow'  control  on  the  conical  forebody  cross-section  using  several 
pulsed  dc  plasma  actuators  show  considerable  improvement  for  separated  neutral  gas 
How  at  17.5  deg  angle  of  attack.  In  general,  two-dimensional  details  are  limited  in 
exploring  discharge  constrictions  and  instability  structures,  and  thus  the  three- 
dimensional  analysis  of  force  generation  mechanism  remains  necessary. 

As  an  extension  of  our  two-dimensional  model  ',  in  this  effort  w  e  cast  the  drift 
diffusion  discharge  physics  into  a  three-dimensional  finite  element  framework.  We  apply 


the  formulation  for  modeling  a  set  of  pulsed  dc  plasma  actuators  over  a  three-dimensional 
conical  forebody  using  plasma  actuator  operating  in  weakly  ionized  oxygen  gas.  The 
section  11  gives  details  of  air  (N2/O2)  chemistry'.  The  third  section  provides  details  of  the 
computational  geometry,  initial  and  boundary'  conditions  and  numerical  details. 
Simulation  of  the  plasma  discharge  physics  with  applied  pulsed  dc  potential  is  described 
in  section  IV.  Conclusions  are  drawn  in  the  final  Section  V. 


II.  AIR  ELECTROCHEMISTRY 

For  the  air  chemistry'  we  neglect  the  metastable  species  along  with  N4  and  O4  due  to 
their  extremely  high  recombination  rates.  Also,  the  numerical  complexity’  is  further 
simplified  by  excluding  nitrous  oxide  at  this  stage.  The  model  equations  governing 


chemistry  of  discharge  are  as  follows.3 

(1)  Nitrogen  model 

(a)  Ionization  e+N2  — >  AC  +2e.  (la) 

(b)  Dissociation  e  +  AC  — >  2N  +  e ,  (lb) 

(c)  Dissociative  recombination  e  +  N£  — »  AC,  (lc) 

(2)  Oxygen  model 

(a)  Ionization  e  +  02  — >  Oj  +2e,  (2a) 

(b)  Dissociation  e  4-  02  — >  2 O  +  e,  (2b) 

(c)  Dissociative  attachment?  4-  Oz  — >  0_  +  O ,  (2c) 

(d)  Recombination  <9"+02+  — >0  +  02<  (2d) 

(e)  Dissociative  recombination  0  +  0~  — >  02  4-  e ,  (2e) 


The  drift-diffusion  form  of  continuity  and  Poisson  equations  for  the  electrons,  ions  and 
neutrals  are  as  follows: 

3n 
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with  momentum  flux  na\a  =-sgn {e)nafiaV$-  DaS/nti,  and  electric  field  E  =  -V^  . 


The  electron  temperature  is  calculated  from  E  =  kBTe/(Vne/ne).  which  is  obtained 
assuming  an  initial  Boltzmann  distribution,  ne  oc  cxp(e<fi /kf]Te).  Various  rate  coefficients 
k i,  A’2, ...^8  related  to  in  Eqs.  (la)  through  (2e)  are  obtained  from  Kossyi  et  al.  2  as 
functions  of  electron  temperature.  The  secondary  electron  emission  from  the  exposed 
surface  is  taken  as  a  function  of  incident  electron  energy  3  and  remains  small,  i.e.,  less 
than  10*  .  No  material  sputtering  of  the  surface  is  considered.  In  Eqns.  (3)-(l  1),  /7f/and  va 
arc  density  and  velocities  of  species  a.  The  mobilities^  and  diffusion  rates  D(I  are  taken 

from  Ellis  et  al  4  Details  of  air  chemistry  simulation  results  considering  the  bulk  density 
of  the  air  as  1.3  kg/m  and  an  atmospheric  ratio  of  3.6  for  nitrogen  to  oxygen  gas 
molecules  arc  described  in  our  publication.  The  three  dimensional  discharge  over 
forcbody  cone  is  simulated  using  oxygen  electrochemistry  where  the  bulk  density  of 
oxygen  is  taken  to  be  0.35  kg/m  at  atmospheric  pressure. 


III.  NUMERICAL  MODEL 


Figure  2  shows  schematic  of  a  pulsed  dc  plasma  actuator.  The  simulated  region  is 
60°  segment  of  a  cone.  The  segment  extents  from  z  =  0.06  cm  to  z  =  4.0  cm.  The  surface 
containing  the  electrodes  and  dielectric  is  given  by  x  -  0.16sin(360-a)  cm  to 
.v  =  0.5sin(360-ar)  cm  and  v -0.16cos(360-a)  cm  to  v  =0.5cos(360-a)  cm.  where  a 
varies  from  0C  to  60°.  The  domain  of  simulation  consists  of  a  2  mm  thick  air  packet 
coating  over  this  curved  surface.  The  domain  size  is  sufficient  to  capture  the  physics  of 
plasma  dynamics. 


To  Ground 


Figure  2.  Schematic  of  three-dimensional  plasma  discharge  on  a  conical  forcbody  using 
pulsed  dc  actuators  operating  in  oxygen  gas. 

We  study  only  plasma  dynamics  in  the  present  work.  Present  domain  size  is  not 
sufficient  to  capture  dynamics  of  air  flow  (for  Navier-Stokes  equations)  which  needs  a 
bigger  domain  size.  The  pulsed  dc  plasma  actuator  consists  of  four  electrodes  separated 
by  a  dielectric.  The  powered  electrodes  are  marked  by  P  and  grounded  electrodes  by  G 
the  Figure  2.  The  electrodes  are  exposed  to  oxygen.  Thickness  of  each  electrode  is  2.0 
mm  along  the  z-direction.  The  powered  and  grounded  electrodes  are  separated  by 
dielectric.  Thickness  of  dielectric  separating  the  electrodes  is  2.0  mm  along  the  z- 
direction.  The  powdered  electrodes  are  from  z  =  1 .6  cm  to  z  =  1 .8  cm  and  z  =  2.4  to  z  =  2.6 
cm.  The  dielectric  surface  between  the  electrodes  is  from  z  =  1.8  cm  to  z  =  2.0  cm,  z  - 
2.2  cm  to  z  =  2.4  cm  and  z  =  2.6  cm  to  z  =  2.8  cm.  The  grounded  electrode  is  from  with  z 
=  2.0  cm  to  z  =  2.2  cm  and  z  =  2.8  to  z  -  3.0  cm.  We  have  chosen  the  relative  dielectric 
constant  er  =  10.  The  domain  is  filled  with  air  of  relative  dielectric  constant  er  -  1.0055. 
The  height  of  the  electrodes  is  assumed  negligible.  A  pulsed  dc  voltage  of  (j)  =  <t>o  sin4 


{2nfi)  volts  is  applied  to  the  powered  electrode.  Such  pulsed  dc  voltage  can  be  produced 
by  commercially  available  function  generators  and  can  be  amplified.  We  have  chosen 
with  <)>()=  1000  volt,  and/=  5  kHz. 

Initial  conditions  are  as  follows:  The  all  initial  particle  concentrations,  except 
those  of  the  electrons,  and  oxygen  molecules,  are  taken  zero.  Initial  atmospheric  oxygen 
molecule  density  is  taken  as  1026/m3,  and  the  electron  density'  is  taken  as  109/m  Initial 
pulsed  dc  potential  is  zero.  These  initial  conditions  have  been  chosen  to  match  with  the 
realistic  atmospheric  conditions. 

The  boundary  conditions  for  the  Poisson’s  equation  are  as  follow's:  The  potential 
difference  with  bo  =  1000  volt  is  applied  to  the  first  electrode  of  the  pair.  The  other 
electrode  of  the  pair  is  grounded.  The  powered  electrodes  are  marked  by  P  and  grounded 
electrodes  by  G  the  Figure  2.  The  powered  and  grounded  electrodes  are  separated  by 
dielectric.  Electric  insulation  condition  (normal  component  of  electric  field  equal  to  zero) 
is  assumed  at  outer  boundaries  of  the  domain  and  at  the  dielectric  surface. 

Boundary  conditions  related  to  oxygen  species  continuity'  equations  are  as 
follows:  the  currents  flow  normal  only  to  the  rf  electrode  since  it  is  an  equipotential 
surface.  Homogeneous  Neumann  conditions  are  applied  to  the  outer  edges  of  the  domain 
and  electric  insulation  is  assumed  at  the  surface  of  the  dielectric.  The  normal  currents  at 
the  surface  of  dielectric  are  taken  zero. 

The  ionized  oxygen  gas  is  numerically  modeled  using  a  self-consistent  three- 
dimensional  Galerkin  variational  formulation  based  finite-element  method  to  obtain 
electron,  ion  and  neutral  species  densities  of  oxygen,  and  the  electric  potential 
distribution.  Finite  element  techniques3,4’14’18  are  well  known  for  their  adaptability  to 
arbitrary  multidimensional  geometries  and  accurate  imposition  of  complicated  boundary 
conditions.  Here,  the  equation  sets  (3):  (11)  can  be  written  with  operator  L  as  L  (U)  =  0 
where  lT  contains  state  variables  like  species  number  densities  and  potential  Multiplying 
with  a  permissible  test  function  rj  and  integrating  over  the  spatially  discretized  domain  Q, 
the  variational  statement  results  in  the  weak  form 


(12) 


for  a  discretization  h  of  Q  =  U^.  and  3C  is  the  non-overlapping  sum  over  the  elements. 

Application  of  Green-Gauss  divergence  theorem  “weakens”  the  order  of  derivatives  in 
(12)  and  yields  natural  homogenous  Neumann  boundary  conditions.  The  surface  integral 
thus  resulting  in  Eqn.  (12)  contains  the  (un)known  boundary  fluxes  wherever  fixed  or 
flux  boundary  conditions  are  enforced  accurately.  Thus,  for  example,  the  Galerkin  form 
of  Eqn.  (11)  becomes, 

/  i  j  \ 

X1  +e  -e  j  t]t]'d:{n0. },  +  e },  -  J  tjs  ~±{</>}e  =  /•„ 

*  c*~  n,  a,  ’  n,  a,  ^,xir  c*z  J 

where  F $  is  the  solution  residual. 

The  Jacobian  matrix  J=[dF/d\i ]  in  the  global  [J].{fU)  =  -{F}  is  resolved  using  LU- 

decomposition  scheme  for  updating  change  in  discretized  solution  vector  U  at  each 
iteration.  The  terminal  non-linear  ordinary  differential  equation  (ODE)  systems  are 
solved  using  implicit  Euler  method  for  temporal  evolution  and  N-R  iterative  algorithm 
for  the  non-linear  matrix  algebra.  The  convergence  criterion  for  all  variables  at  any 
iteration  is  set  at  10~4.  Solution  stability  is  ensured  by  appropriate  selection  of  adaptive 
time  marching  step  size  (with  the  smallest  step  of  one  hundredth  of  the  dielectric 
relaxation  timescale)  and  the  introduction  of  artificial  diffusion.  The  number  of  degrees 
of  freedom  employed  is  123816. 

Independent  of  the  physical  dimension  of  the  working  domain  12,  and  for  general  forms 
of  the  flux  vectors,  the  semi-discretized  weak  statement  always  yields  an  ordinary 
differential  equation  (ODE)  system  that  is  fully  discretised  using  a  ^-implicit  or  r-step 
Runge-Kutta  type  time  integration.  The  terminal  ODE  is  usually  solved  using  a  Newton- 
Raphson  scheme  for  U(/)- 

U'/J,  =  U'r+I  +  AU'  =  Ur  +XU',+1-  where 

P-0 

AU'  =-[M+<9A/(dR/dU)]'’R(U) 


(13) 


Here,  a  <9— implicit  time  marching  procedure  is  employed.  In  (13).  M  -  Se(Me)  is  the 
“mass”  matrix  associated  with  element  level  interpolation,  R  carries  the  element 
convection,  diffusion  and  source  information.  The  calculation  of  the  "Jacobian”  SR/SU 
and  inversion  of  the  M+#A/(3R/dU)  matrix  with  sufficient  accuracy  is  obviously  a 
numerical  challenge.  However,  unlike  the  traditional  finite  difference/volume  methods, 
the  present  FE  algorithm  allows  one  to  simulate  the  system  simultaneously  without 
requiring  any  sub-iteration  for  the  Poisson  solver. 

In  order  to  resolve  the  plasma  details  in  three-dimensions  it  is  crucial  to  resolve 
wall  sheath  region.  So,  we  used  the  smallest  element  thickness  of  the  order  of  Debye 
length  around  the  powered  electrode.  The  mesh  density  was  proportional  to  the  electric 
field  intensity  in  all  three-directions.  Note  in  this  case  we  are  not  attempting  to  calculate 
fluid  flow  which  would  require  a  very  large  domain  with  sub-mm  scale  flow  mesh  to  be 
overlapped  on  the  micron  scale  plasma  mesh.  In  such  case,  in  three-dimensions  several 
hundreds  of  millions  of  degrees  of  freedoms  will  be  solved  mandating  some  multi  level 
parallel  coding  approach  and  weeks  of  simulation  time  to  cover  all  time  scales.  The  level 
of  computational  challenge  will  even  increase  further  as  we  include  more  complexity'  of 
finite  rate  air  chemistry  with  large  species  array.  We  contain  such  challenges  by  solving 
drift-diffusion  approximated  plasmas  incorporating  only  electron  and  four  species  of 
oxygen  gas.  The  select  set  of  equations  from  Eqn.  (1):  (11)  are  used  for  modeling  oxygen 
chemistry  as  given  in  our  recent  paper.  3 


IV.  RESULTS  AND  DISCUSSION 

Discharge  starts  in  the  oxygen  as  we  switch  on  the  pulsed  dc  voltage.  Different 
ion  and  neutral  species  are  formed  through  ionization,  dissociation  and  dissociative 
attachment  etc.  Recombination  also  occurs.  Full  chemistry  of  different  species  formation 
and  recombination  is  described  by  Equation  (2).  The  simulation  results  at  the  peak  of  the 
10th  cycle  of  the  pulsed  dc  are  given  in  Figures  3-7. 

Figure  3(a)  and  3(b)  showr  electric  potential  and  electric  field  components  as  a 
function  of  z  at  a  =  15°  and  45°  for  a  pulsed  dc  plasma  actuator.  As  mentioned  earlier 
that  very  first  electrode  from  the  left  is  powered  electrode.  Then  there  is  dielectric  of  2 


mm  thickness  and  then  a  grounded  electrode.  The  sequence  is  repeated  for  second  set  of 
electrodes  and  dielectric.  The  variation  of  various  components  of  electric  field  on  the 
surface  follows  from  applied  of  electric  potential.  Electric  field  is  given  byE  =  -V^. 
Since  the  gradient  of  electric  field  is  highest  at  the  edges  the  electrodes,  the  value  of 
electric  field  is  highest  there.  The  potential  does  not  change  with  z  for  the  powered  and 
grounded  electrodes.  The  gradient  of  electric  potential  with  respect  to  z  is  zero.  The  value 
of  Ez  is  zero  at  the  electrodes  surfaces.  The  potential  is  not  constant  with  respect  to  x  or  y, 
hence  Ex  and  Ev  are  not  zero  at  the  electrode  surfaces.  The  curvature  of  the  electrodes 
with  lower  values  of  z  is  high;  the  electric  field  is  high  for  these  electrodes.  Higher 
ionization  occurs  for  these  electrodes. 

Figure  4  (a)  shows  spatial  distribution  of  electron  density  for  pulsed  dc  plasma 
actuator.  Chemistry  of  electrons  formation  is  given  by  set  of  Equation  (2a).  Spatial  and 
temporal  density  profile  of  electrons  is  governed  by  continuity  equations  coupled  with 
Poission  equation.  Electron  density  rises  up  to  a  level  4  x1014/m  .  Electrons  respond  to 
electric  field  very  fast  due  to  low  mass  and  high  mobility.  The  electrons  are  repelled  from 
the  area  of  grounded  electrode  and  are  attracted  towards  powered  electrode  where  they 
are  lost.  Figure  5  (a)  shows  slices  of  density  of  electrons  at  the  positive  peak  of  the  pulsed 
dc  of  tenth  cycle.  The  density  decreases  sharply  with  the  increase  in  the  distance  from  the 
surface  of  the  electrodes  and  dielectric. 


Z  (cm)  Z  (cm) 


(a)  a  =  15"  (b)  a  =  45' 

Figure  3.  Electric  potential  and  electric  field  components  as  a  function  of  z. 

The  oxygen  ions  0“  are  produced  through  dissociative  attachment  of  oxygen 
molecules  with  electrons  governed  by  the  Eqn.  (2c).  Its  density  is  governed  by  continuity 


Eqn.  (6)  coupled  with  other  equations.  The  density  of  oxygen  ions  O  grows  to  a  level  of 
nearly  3  ><10I5/m3.  Figure  4  (b)  shows  density  profile  of  negative  oxygen  ions  O  at 
positive  peak  of  the  pulsed  dc  of  tenth  cycle.  The  negative  oxygen  ions  play  an  important 
role  in  oxygen  chemistry.  The  oxygen  ions  O”  are  repelled  from  the  area  of  grounded 
electrode  and  are  attracted  towards  powered  electrode  under  the  influence  of  pulsed  dc. 
At  the  electrode  oxygen  ions  O”  become  a  part  of  normal  current.  The  oxygen  ions 
O  are  much  heavier  than  that  of  electrons  and  its  mobility  is  low;  response  is  slow. 

The  edge  of  the  powered  electrode  is  at  z=1.8  cm.  The  peak  of  density  is  close  to 
the  edge  of  the  powered  electrode.  Electric  field  is  highest  at  the  edges  hence  higher 
ionization  occurs  at  these  locations.  The  curvature  of  the  electrodes  at  lower  values  of  z  is 
higher.  The  electric  field  is  higher  for  higher  curvature.  This  increases  ionization  for  the 
electrodes  with  lower  values  of  z.  Figure  5  (b)  shows  slices  of  density  of  negative  oxygen 
ions  0“  at  positive  peak  of  the  pulsed  dc  of  tenth  cycle  at  electrode  edges.  The  density 
decreases  sharply  with  the  increase  in  the  distance  from  the  edge  of  the  electrodes  or 
dielectric. 

Response  of  negative  species  to  the  electric  field  is  opposite  to  that  of  positive 
ions.  Ion  mass  and  motilities  are  also  different  for  different  oxygen  species  and  electrons, 
hence,  their  velocities  are  different.  This  gives  rise  to  charge  separation  given  by 
n  =n  -n  —n  .  Figure  4  (d)  shows  density'  profile  of  charge  separation  at  positive 

peak  of  the  tenth  cycle.  It  resembles  that  of  density  profiles  of  positive  oxygen  ions.  The 
densities  of  charge  species  (such  as  electrons  and  02  )  are  small  in  comparison  to  O" 
ions.  The  value  of  charge  separation  q  is  nearly  equal  to  density  of  (X  ions.  Figure  5  (d) 
shows  slices  of  density  of  charge  separation  at  the  edges  of  electrodes  at  positive  peak  of 
the  pulsed  dc  of  tenth  cycle.  The  peak  of  charge  separation  is  very  close  to  the  common 
edge  of  the  powered  electrode  and  dielectric.  Strength  and  direction  of  electric  field  and 
charge  of  various  species  is  responsible  for  this  location  of  densities.  Charge  separation 
gives  rise  to  a  self-generated  electric  field  through  Poisson’s  Equation.  The  density 
decreases  sharply  with  the  increase  in  the  distance  from  the  edges  of  the  electrodes  or 


dielectric.  The  localized  hotspots  near  the  electrodes  are  similar  to  that  seen  in 
experiments  and  may  be  signatures  of  DC  corona  instability. 
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Figure  4.  Surface  plot  of  density  (nT  )  of  (a)  electrons,  (b)  negative  oxygen  ions  0~  , 
(c)  positive  oxygen  ions  O*  and  (d)  charge  separation  nq  -  n  y  -  itc  -  n  . 


Figure  5.  Slices  of  density  (m3)  of  (a)  electrons,  (b)  negative  oxygen  ions  O  ,  (c) 
positive  oxygen  ions  Of  and  (d)  charge  separation  nq  =n  -nc  —  n  . 


Figures  6  (a)  and  6  (b)  show  density  of  various  charge  speeies  at  a  =  15°  and  45° 
as  a  function  of  z,  respectively.  These  figures  give  us  numerical  estimation  of  variation  of 
density  of  different  charge  speeies  along  z.  Since  the  discharge  is  not  uniform  along  the 
edges  of  the  electrons,  variation  of  charge  densities  in  Figs.  6  (a)  and  5(b)  are  not  same. 

xIO15  xIO15 


z  (cm) 

(a) 

xIO16  xIO15 


z  (cm) 

(b) 

Figure  6.  Density  of  various  charge  speeies  at  (a)  a  =  15"  and  (b)  a  -  45". 

Figure  7  shows  spatial  distribution  of  vectors  of  time  average  of  eleetrodynamic 
force  F  =  e(w  -nr-«0  )E  for  pulsed  dc  plasma  actuator.  The  force  is  strong  elose  to  the 

common  edges  of  dielectric  and  electrodes  due  to  high  eleetrie  field  and  ionization.  The 


direction  of  time  averaged  force  is  in  positive  z-direction  from  powered  electrode  to 
grounded  electrode  (from  anode  to  cathode).  This  may  be  useful  in  separation  control. 
The  solution  will  be  sensitive  to  the  initial  conditions;  we  have  taken  only  one  initial 
condition  which  is  close  to  the  atmospheric  condition  We  have  investigated  the  role  of 

negative  oxygen  ions  02  in  air  chemistry  and  found  that  the  density  of  negative  oxygen 

ions  02  remains  lower  than  that  of  positive  oxygen  ions  02 .  The  positive  oxygen  ions 

02  play  more  important  role  compared  to  02  .  These  results  have  been  described  in  our 
recent  study  3<\ 


Figure  7.  Vector  plot  of  electrodynamic  force  F  =  e(nQ.  -  nc  —  it  )E  in  perspective 
shows  three-dimensional  nature  of  the  force  distribution. 

V.  CONCLUSIONS 

For  the  first  time,  the  three-dimensional  discharge  about  a  plasma  actuator 
operating  in  oxygen  under  pulsed  dc  power  has  been  studied  using  a  self-consistent 
multibody  modeling  of  oxygen,  plasma  and  dielectric.  Chemistry  of  formation  of 
different  species  of  oxygen  has  been  taken  into  account.  Species  with  very  high 
recombination  rates  and  nitrous  oxide  have  been  neglected  for  simplicity.  Continuity 
equations  governing  densities  of  electrons,  ions  and  neutral  species  of  nitrogen  and 
oxygen  are  solved  with  Poisson's  equation  using  a  three-dimensional  finite  element  based 
formulation  of  plasma  to  obtain  spatial  and  temporal  profiles  of  densities  of  species,  and 
electric  field.  The  electric  field  and  ionization  is  the  highest  close  to  the  common  edges  of 
electrodes  and  dielectric.  Oppositely  charged  species  move  in  opposite  directions  due  to 
applied  pulsed  dc  potential  which  gives  rise  to  charge  separation.  Poisson’s  equation 
governs  electric  field  generated  due  to  charge  separation.  The  density  of  separated  charge 


n  =  (//  -nc -n  ) ,  and  the  electrodynamic  force  per  volume  F  =  e(n().  -nc -nQ.  )E 

have  been  obtained.  Numerical  predictions  show  DC  corona  type  plasma  instability 
similar  to  experimental  observations.  The  time  average  of  the  force  shows  mostly 
acceleration  in  the  forward  direction  (from  anode  to  cathode)  above  the  actuator.  This 
results  in  a  moving  wave  of  plasma  over  the  surface  in  forward  direction  which  can  find 
application  in  flow  control.  A  pulsed  dc  plasma  actuator  may  be  used  as  an  alternative  to 
the  dielectric  barrier  discharge  plasma  actuator. 

As  a  next  step  we  would  like  to  validate  our  numerical  model  integrating  with 
diagnostic  experiments  as  follows.  (1)  Study  the  barrier  discharge  plasma  actuators  in 
three-dimensions  using  first  principles  models;  (2)  Utilize  the  optical  and  spectroscopic 
diagnostics  to  improve  the  air  chemistry  models;  (3)  Introduce  design  modifications  for 
improved  actuator  performance;  and  (4)  Transition  developed  formulations  to  the 
RBAC/AFRL  researchers  in  a  reduced  order  form  for  incorporation  into  their  fluid 
dynamic  codes  in  a  loosely  coupled  manner.  The  research  will  culminate  into  the  most 
sophisticated  modeling  and  the  most  detailed  characterization  of  plasma  based  flow 
control  devices  available  to-date.  Gaining  insight  of  plasma  actuators  will  help  robust, 
scalable  actuation  of  transitional  and  turbulent  flows. 

Acknowledgment/Disclaimer 

This  work  was  sponsored  by  the  Air  Force  Office  of  Scientific  Research,  USAF,  under 
grant  number  FA9550-07-1 -01 3 1  monitored  by  Dr.  John  Schmisseur  and  LtCol  Rhett 
Jefferies.  The  views  and  conclusions  contained  herein  are  those  of  the  authors  and 
should  not  be  interpreted  as  necessarily  representing  the  official  policies  or  endorsements, 
either  expressed  or  implied,  of  the  Air  Force  Office  of  Scientific  Research  or  the  U.S. 
Government.  The  authors  acknowledge  thoughtful  discussions  with  Drs.  Datta  Gaitondc 
and  Miguel  Visbal. 


REFERENCES 

1  J.R.  Roth.  D.M.  Sherman.  S.R.  Wilkinson.  AIAA  J.  38,  7  (2000). 

2  C.  L  Enloe,  T.  E.  McLaughlin,  G.  I.  Font,  and  J.  W.  Baughn,  44th  AIAA  Aerospace 
Sciences  Meeting  and  Exhibit,  Reno,  Nevada,  9-12  January  2006  (AIAA,  Washington, 
D.C.  2002),  Paper  No.  AIAA  2006-166. 

S.Roy  and  D.V.Gaitonde.  43rd  AIAA  Aeropsace  Sciences  Meeting,  Reno  NV,  10-13 
January,  2005  (AIAA,  Washington,  D.C.  2002),  Paper  No.  AIAA-2005-0160. 

4  S.Roy,  K.P.  Singh  and  D.V.Gaitonde,  Appl.  Phys.  Lett.  88,  121501  (2006). 

J.  Poggie.  36th  AIAA  Plasmadynamics  and  Lasers  Conference,  Toronto.  Ontario. 
Canada.  6-  9  June  2005  (AIAA.  Washington.  D.C.  2002),  Paper  No.  AIAA  2005-5303. 
"K.P.  Singh  and  S.  Roy,  Modeling  plasma  actuators  with  air  chemistry  for  effective  flow 
control.  J.  App.  Phys.  101.  123308  (2007). 

B.  Jayaraman,  S.  Thakur.  W.  Shyy,  J.  Heat  Transfer  129.  517  (2007). 

*  A.V.  Likhanskii.  M.N.  Shneider,  S.O.  Macheret  and  R.B.  Miles,  Phys.  Plasmas  14. 
073501  (2007). 

9  T.C.  Corke,  M.L.  Post  and  D.M.  Orlov,  Prog.  Aerospace  Sci.  43,  193  (2007). 

10  G.  1.  Font,  C.  L.  Enloe,  T.E.  McLaughlin,  and  D.  Orlov,  45th  AIAA  Aerospace 
Sciences  Meeting  and  Exhibit  Reno,  Nevada  2007  (AIAA,  Washington,  D.C.  2002), 
Paper  No.  AIAA  2007-188. 

J.  W.  Gregory,  C.  L.  Enloe,  G.  I.  Font,  and  T.  E.  McLaughlin,  45th  AIAA  Aerospace 
Sciences  Meeting  and  Exhibit  Reno,  Nevada  2007  (AIAA,  Washington,  D.C.  2002). 
Paper  No.  AIAA  2007-185. 

"  S.  Pavon,  J-L  Dorier,  Ch  Hollenstein,  P.  Ott  and  P.  Leyland,  J.  Phys.  D:  Appl.  Phys. 
40,  1733  (2007). 

13  C.  O.  Porter,  J.  W.  Baughn,  T.  E.  McLaughlin,  C.  L.  Enloe,  and  G.  1.  Font,  44th  AIAA 
Aerospace  Sciences  Meeting  and  Exhibit  Reno,  Nevada,  9-12  January'  2006  (AIAA. 
Washington,  D.C.  2002),  Paper  No.  AIAA  2006-104. 

14  K.P.  Singh  and  S.  Roy,  J.  Appl.  Phys.  103,  013305  (2008). 

15  M.P.  Patel.  T.T.  Ng  and  S.  Vasudevan,  T.C.  Corke,  M.L.  Post,  T.E.  McLaughlin  and 
C.F.  Suchomel,  45th  AIAA  Aerospace  Sciences  Meeting  and  Exhibit.  Reno,  NV,  8-11 
January  2007  (AIAA.  Washington,  D.C.  2002),  Paper  No.  AIAA  2007-635. 


16  T.  C.  Corke,  B.  Mertz,  M.  P.  Patel,  44th  A1AA  Aerospace  Sciences  Meeting  and 
Exhibit  Reno,  Nevada,  9-12  January  2006  (AIAA,  Washington,  D.C.  2002),  Paper  No. 
A1AA  2006-1208. 

J.  R.  Roth  and  X.  Dai,  44th  AIAA  Aerospace  Sciences  Meeting  and  Exhibit,  Reno, 
Nevada,  9-12  January  2006  (AIAA,  Washington,  D.C.  2002),  Paper  No.  AIAA  2006- 
1203. 

*  K.P.  Singh  and  S.  Roy,  Appl.  Phys.  Lett.  91,  081504  (2007). 

19  C.  Baird,  C.  L.  Enloe,  T.  E.  McLaughlin,  and  J.  W.  Baughn,  Proceedings  of  the  43rd 
Aerospace  Sciences  Meeting,  Reno,  Nevada,  2005  (AIAA,  Washington,  D.C.,  2002), 
Paper  No.  AIAA-2005-0565. 

20  W.  Shyy,  B.  Jayaraman,  and  A.  Anderson,  J.  Appl.  Phys.  92,  6434  (2002). 

21  S.  Roy,  Appl.  Phys.  Lett.  86,  101502  (2005).  K.  P.  Singh  and  S.  Roy,  J.  Appl.  Phys.  98, 
083303  (2005):  B.  Jayaraman,  S.  Thakur,  and  W.  Shyy,  J.  Heat  Transfer  129.  517  (2007); 
A.  V.  Likhanskii,  M.  N.  Shneider,  S.  O.  Macheret,  and  R.  B.  Miles,  Phys.  Plasmas  14, 
073501  (2007). 

22  D  Gaitonde,  M.  Visbal.  and  S.  Roy,  ASME  Joint  U.S. -European  Fluids  Engineering 
Summer  Meeting,  Miami,  2006  (AIAA,  Washington,  D.C.,  2002),  Keynote  Lecture 
FEDSM2006-98553;  M.  Visbal,  D.  Gaitonde,  and  S.  Roy,  Proceedings  of  the  Fluid 
Dynamics  and  Flow  Control  Conference,  San  Francisco.  June  2006  (AIAA,  Washington, 
D.C.,  2002),  Paper  No.  AIAA-2006-3230. 

23  K.  P.  Singh,  S.  Roy,  and  D.  Gaitonde,  Plasma  Sources  Sci.  Technol.  15,  735  (2006). 
*4K.P.  Singh  and  S.  Roy,  Applied  Physics  Letters,  91, 081504  (2007). 

O.A.S.  Kandil,  H.  Hazem  and  C.H.  Liu,  AIAA  Atmospheric  Flight  Mechanics 
Conference,  Hilton  Head  Island,  (SC,  Washington,  1992),  AIAA  1992-244. 

26  B.  Jayaraman.  W.  Shyy,  Prog.  Aerospace  Sci.  (in  press,  2007). 

27  Celik,  Zcki  Z.  Roberts,  Leonard,  AIAA- 1992-4430  AIAA  Atmospheric  Flight 
Mechanics  Conference,  Hilton  Head  Island,  (SC,  Washington,  1992),  Paper  No.  AIAA 
1992-293. 

28  N.  Pedreiro.  S.M.  Rock,  Z.Z.  Celik  and  L.  Roberts,  Journal  of  Aircraft  35,  69  (1998). 


29  V.  Shalaev,  A.  Fedorov,  N.  Malmuth.  V.  Zharov,  and  Ivan  Shalaev,  41st  Aerospace 
Sciences  Meeting  and  Exhibit,  Reno.  Nevada,  6-9  January  2003,  (A1AA.  Washington. 
D.C.  2002).  Paper  No.  AIAA  2003-34. 

30  A.  A.  Maslov.  B.  Y.  Zanin,  A.  A.  Sidorenko,  B.  V.  Postnikov.  V.  P.  Fomichcv,  A.  D. 
Budovsky,  and  N.  Malmuth,  43r1  AIAA  Aerospace  Sciences  Meeting  and  Exhibit.  Reno, 
Nevada,  2005.  (AIAA,  Washington.  D.C.  2002),  Paper  No.  AIAA-2005-400. 

31  K.P.  Singh  and  S.  Roy,  J.  App.  Phys.  101,  093301  (2007). 

2 1.  A.  Kossyi,  A.  Yu  Kostinsky,  A. A.  Matveyev  and  V.P.  Silakov,  Plasma  Sources  Sci. 
Techol.  1.  207  (1992). 

1  S.  Roy  and  B.P.  Pandey,  Phys.  Plasmas,  9,  4052  (2002). 

4  H.W.  Ellis,  R.  Y.  Pai.  E.  W.  McDaniel,  E.  A.  Mason,  and  L.  A.  Viehland,  At.  Data 
Nucl.  Data  Tables,  17.  177  (1976). 

35  K.P.  Singh  and  S.  Roy,  Journal  of  Applied  Physics,  103  (10)  103303  (2008). 

36  K.P.  Singh  and  S.  Roy.  Applied  Physics  Letters,  92  ( 1 1 )  111  502  (2008). 


Personnel  Supported  During  Duration  of  Grant 


Richard  Anderson 
Chin-Cheng  Wang 
Kunwar  P.  Singh 
Sub  rata  Roy 


Graduate  Student,  University  of  Florida,  Gainesville 
Graduate  Student,  University  of  Florida  Gainesville 
Post-doctoral  Associate,  University  of  Florida,  Gainesville 
Associate  Professor,  University  of  Florida,  Gainesville 


Relevant  Publications  (2007-2008) 

Wang,  C.C.  and  S.  Roy,  “Electrodynamic  enhancement  of  film  cooling  of  turbine 
blades,"  J.  Appl,  Phys.  104  (7)  Scheduled  for  Oct  2008. 

Roy,  S.  and  C.C.  Wang,  2008,  “Plasma  actuated  heat  transfer,”  Applied  Physics  Letters 
92,  231501. 


Singh  K.P.  and  S.  Roy,  2008,  “Force  generation  due  to  three  dimensional  plasma 
discharge  on  a  conical  forebody  using  pulsed  dc  actuators,”  Journal  of  Applied  Physics, 
103  (10)  103303. 

Singh  K.P.  and  S.  Roy,  2008,  “Physics  of  plasma  actuator  operating  in  atmospheric  air," 
Applied  Physics  Letters,  92  (1  1)  11  1502. 


Singh  K.P.  and  S.  Roy,  2008,  “Force  approximation  for  a  plasma  actuator  operating  in 
atmospheric  air,"  Journal  of  Applied  Physics,  103  (1)  013305. 

Singh  K.P.  and  S.  Roy,  2008,  “Flow  control  on  a  three  dimensional  forebody  using 
pulsed  dc  plasma  actuator  operating  in  air,”  AIAA-2008-1 102,  46th  AIAA  Aerospace 
Sciences  Meeting  and  Exhibit 

Singh  K.P.  and  S.  Roy,  2007,  “Impedance  matching  for  an  asymmetric  dielectric  barrier 
plasma  actuators,”  Applied  Physics  Letters,  91  (8)  081504. 

Singh  K.P.  and  S.  Roy,  2007,  “Modeling  plasma  actuators  with  air  chemistry  for  effective 
flow  control,”  Journal  of  Applied  Physics,  101  (12)  123308. 

Singh  K.P.  and  S.  Roy,  2007,  “Vortical  flow  control  on  a  conical  forebody  cross-section 
using  an  array  of  pulsed  dc  actuators,”  Journal  of  Applied  Physics,  101  (9)  093301. 

Roy,  S.,  K.P.  Singh,  and  D.  Gaitonde,  2007,  “Air  plasma  actuators  for  effective  flow 
control,”  AIAA-2007-1 84,  45th  Aerospace  Sciences  Meetings,  Reno,  NV. 


